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Abstract. We introduce EXOFIT, a Bayesian tool for estimating orbital 
parameters of extra-solar planets from radial velocity measurements. EXOFIT 
can search for either one or two planets at present. EXOFIT employs Markov 
Chain Monte Carlo method implemented in an object oriented manner. As 
an example we re-analyze the orbital solution of IID155358 and the results are 
compared with that of the published orbital parameters. In order to check the 
agreement of the EXOFIT orbital parameters with the published ones we exam- 
ined radial velocity data of 30 stars taken randomly from www. exoplanet. eu. We 
show that while orbital periods agree in both methods, EXOFIT prefers lower 
eccentricity solutions for planets with higher (e > 0.5) orbital eccentricities. 



Introduction 



More than a decade of extensive search for extra-solar planets has resulted in 
nearly 300 planets. Majority of the contribution to the extra-solar planet count 
comes from radial velocity m ethod. Radial velocity dat a is traditionally ana- 
lyzed first by a periodogram ( Lomb I [19731 : IScargldfl982l ) for the orbital period 
and then for the other orbital parameters using conventional optimization meth- 
ods. Bayesian methods fo r the esti r nation of o rbital parameters of extra-solar 
planets were introduced by [Gregory (|2005al ) and lFordI (j2005l ). Their work shows 
that these methods provide us a framework to tackle the problems associated 
w ith the traditional met hods in a transparent and robust manner. exofitQ 
of Balan Sz Lahav ( 20081 ) is the first publicly available package for Bayesian esti- 
mation orbital parameters from radial velocity measurements. In this article we 
discuss the application of EXOFIT to the radial velocity data of HD155358 an d 
compare our results with the published orbital solution ( Cochran et al. II2007I ). 
We also show, by analyzing a randomly selected sample of radial velocity data 
from www.exoplanet.eu that EXOFIT disagrees with the published orbital ec- 
centricities in many instances while the orbital periods from EXOFIT matches 
closely with the published ones. 



^ www . star . ucl . ac . uk/ ~ lahav/ exofit .html 
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The rest of the article is organized as follows. In Section [X] we give a 
brief introduction to radial velocity modeling. Section [33 describes Bayesian 
approach to parameter estimation. EXOFIT is introduced in Section \^ and 
its application to the radial velocity data is discussed in section [53 Section [6] 
explains the analysis of radial velocity data of a sample of stars taken form 
www.exoplanet.eu using EXOFIT. We conclude this article in Section [TD and 
provide an outline on planned work. 



2. Modeling of Radial Velocity Data 

Radial velocity data consists of a set of measured radial velocity entries, cor- 
responding time of observation and the uncertainty in each measurem ent. Ob- 
served radial velocity data is modeled by the equation ( Gregorvll2005al ) 



di = i^i + €i + 6, (1) 

where di is the measured radial velocity data for the ith instant of time tj, i^i 
is true radial velocity of the star, ej is the measurement error assigned by the 
observer and 5 represents any unknown noise present (e.g., signal from another 
planet) in the data. For a statistician, (5 is a nuisance parameter. The true 
radial velocity can be simulated using a mathematical model. Disregarding any 
interactions between planets, the radial velocity of a star for a typical n-planet 
model can be approximately written as a linear combination of n single planet 
radial velocities. Thus, 

n 

V = V Ki{sm{fi + Wi) + ei sintUj), (2) 
1=1 

where V, K, f, e, w represent the systematic velocity, amplitude, true anomaly, 
eccentricity and the longitude of periastron respectively. For the full formalism 
see the User's Guide to EXOFIT. 



3. Bayesian Parameter Estimation 

Bayesian paradigm has its origins in an article published posthumously by Rev 
Thomas Bayes in 1763. Since then the theorem has played a central part in 
proba bilisti c inference. For latest examples in cosmology see iFeroz &: Hobson" " 



(l2008l 1 andE ewis Sz Bridle! (j2002l ). Bayesian methods for the estimation o f or- 



bital param eters of extra-solar planets were introduced by [Gregory (|2005al l and 



FordI (|2005[ ) and their research show that this approach has an edge over the 



traditional methods when dealing with for e.g., highly eccentric orbits. These 
methods also provide a straight forward way of dealing with nuisance parame- 
ters and a robust way of estimating uncertainties associated with the estimates 
of orbital parameters. 

Bayes' theorem for a set of parameters in model H and data D can be 
written as, 

, Pr(D|0,i7)Pr(0im 
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In the above equation Pr(0|D, H) is the posterior probabihty distribution of pa- 
rameters, Pr(D|0, H) is the hkehhood of the data, Pr{&\H) represents the prior 
probabihty distribution of the parameters and Pr(D|ff) is caUed the Bayesian 
evidence. For parameter estimation problems we could simply write the above 
equation as 

Pr(eiD, H) oc Pr(D[0, H) Pr(0|i?) (4) 

Computing the right hand side of Equation H] is the central point of any 
Bayesian parameter estimation . Analytical solutions can be derived for some 
special cases. However, in general we use numerical methods to compute the 
posterior distribution. Although many approximation methods exist, this area 
is dominated by Markov Chain Monte Carlo (MCMC) and other sampling meth- 
ods. These rely on the a random walk through the parameter space and make 
use of the fact that posterior density is proportional to the number of points vis- 
ited in the volume considered. We calculate the marginal posterior distribution 
of each parameter by simple plotting a histogram of the final set of samples. 

Bayesian modeling of the problem consists of defining each of the compo- 
nents mentioned in Equation [H This will be discussed from the context of esti- 
mation of orbital parameters when we apply EXOFIT to the published radial 
velocity data of HD155358 in Section [5l 

4. EXOFIT 

EXOFIT is an easy to use documented software for estimating the orbital pa- 
rameters from radial velocity measurements and it is freely available. It is based 
on Bayesian MCMC method and is implemented in an object oriented frame- 
work in CH — h. It can be easily extended to analyze the radial velocity data of 
more than two planets as well as data from transit photometry. Output of EX- 
OFIT is a set of posterior samples. These can be analyzed using any standard 
statistical software. This explained in User's Guide which can be download from 
EXOFIT websit^l. We also provide sample script for R statistical environmen10 
to analyze posterior samples. Improved parameterization for the problem and 
novel sampling techniques for EXOFIT are under development. 

5. Application to Radial Velocity Data of HD155358 

In this section we develop a Bayesian model for the parameter estimation and ap- 
ply EXOFIT to extract the orbital parameters of the companions of IID155358. 
We start by defining the likelihood of the data. As mentioned in Section [23 we 
assume that the data consists of true radial velocity and some noise. The com- 
ponent of noise arising from the known measurement errors is assumed to be 
normally distributed with standard deviation ai for the ith entry in the data. 
The probability distribution for the unaccounted noise component 6 is chosen 
to a Gaussian distribution with finite variance s^. Therefore, the distribution 
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of the combination €i -\- 5 can be considered as a Gaussian with a variance of 



Assuming each measurement error to be independent and since they 
fohow a Gaus sian distribution , the hkeUhood of data can be written as a product 
GaussiansI Gregory 2005al ) where is the number of entries in the data. 
Therefore, 

Fi(D\@,H) =Aexp ^ ^ '■ 



Vi 



i 



where 



(2vr) 



'N/2 



N 



n i'^' + 



-1/2 



(5) 



(6) 



Thus, our parameter space is {V,Ti, Ki,ei,wi,xijT2, K2,e2,W2,X2, s}, first 11 
from the mathematical model and last one representing the nuisance parameter 
6. The parameter is x is defined for computational purposes and marks the 
periastron passage time as function period T. For more details please consult 

EXOFIT User's Guide. 

Table H taken from iBalan fc Lahavl (120081 ) gives the prior probability dis- 
tributions of each parameter in the model. These priors are chosen in such a 
way that they allow likelihood term in the Equation U] to dominate the posterior 
distribution and thus ensuring inference to be drawn from the observed data. 



Table 1. The assumed prior distribution of orbital parameters and their 
boundaries for a 2-planet model. 



Para. 

V{ms-'-) 
Tiidays) 



Prior 



Mathematical Form Min 



Max 



"T" 



ei 

XI 
Tiidays) 



£2 
X2 



Uniform 
Jeffreys 



Ki{ms~^) Mod. Jeffreys 



Uniform 
Uniform 
Uniform 
Jeffreys 



K2{ms ^) Mod. Jeffreys 



Uniform 
Uniform 
Uniform 

Mod. Jeffreys 



^1 (?T^) 



In 



V I<1 

1 
1 

27r 
1 
1 



T2 In I 



In 



' -K20 + -K2max 

1 
1 

2ir 
1 



-2000 
0.2 



0.0 





0.2 



0.0 









2000 
15000 



2000 
1 

27r 
1 

15000 



2000 
1 

2tt 
1 

2000 



Application of E XOFIT revealed the posterior distribution of planets as 
given in Figure 1(a) and the corresponding radial velocity cu rve is shown in 
Figur e 1(b) We compare our results to the published results by ( Cochran et al. I 
20071 ^ in Tabled Although our results looks similar, we notice a noise factor(s) 
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of 5.43 ms~ which indicates the presence of an additional signal in the data. 
We did look for a third planet in the system, but our results were inconclusive. 
This issue can be settled if we have more observations down the line. 




(a) Marginal posterior destribu- (b) Radial velocity curve 

tions 

Figure 1. Marginal postrior distribuotions of orbital orbital parameters of 
the companions to HD155358 are shown in the left panel. Right panel shows 
the corresponding radial velocity curve for HD155358. 



6. EXOFIT vs exoplanet.eu 

In this section we compare the published orbital periods and eccentricities of 30 
extra-solar planets taken randomly from www.exoplanet.eu to that obtained by 
analyzing corresponding radial velocity data with EXOFIT . The results are 
summarized in Figure [2j It can be observed that the orbital periods extracted 
with EXOFIT matches closely with published ones. However, the eccentricities 
show apparent variation from the published values. In fact EXOFIT tends to 
obtain lower eccentricity solutions. This is clearly noticeable for planets with 
eccentricities greater that 0.5 (i.e. ~ 10% of the planet population). Median was 
used as the point estimator for marginal posterior distribution orbital parameters 
of exo-planets extracted with EXOFIT while comparing with the published 
results. This suggests that the orbital eccentricity is poorly constrained in many 
occasions because of the sparse sampling of the data points. 



7. Summary and Future Work 

W e have introduced EXO FIT by reanalyzing the orbital solution of IID155358 
bv ICochran et a~ ( 20071 ). EXOFIT provides a full Bayesian analysis of the 
problem and spits out the marginal posterior distribution of orbital parame- 
ters. We have also compared the orbital parameters obtained by EXOFIT to 
the published orbital solutions of 30 extra-solar planets taken randomly from 
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Table 2. A table showing the summary of the posterior distribution of or- 
bital parameters e xtracted with EXOFI T 2-planet model and the published 
orbital solution bv lCochran et al.1 (|2007( ). Columns 2, 3 and 4 show the pos- 
terior mean (and standard deviation), median (and 25% and 75% quantiles) 
and the maximum a posteriori, i.e. posterior mode (and 68.3% highest density 
regions) respectively. Mass of the star was assumed to be 0.87M©. 



Parameters 


EXOFIT (Mean) 


EXOFIT(Median) 


EXOFIT (Mode) 


Cochran et al. 


V{ms-^) 


11.09 ± 1.37 




11.23l-i^ 




Tiidays) 


195.22 ±0.91 


195.20lg:«^ 


195 1 7+0-87 


195 ± 1.1 




33.56 ± 2.22 


33.59llj« 


00 r:o + 2.07 
OO.JO_2 Q2 


34.6 ± 3 


ei 


0.13 ±0.04 




0.131H! 


0.112 ± 0.037 


■wi (degrees) 


160.78 ± 19.73 


160.46lj2;77 


158.20lf,;^° 


162 ± 20 


XI 


0.38 ±0.07 


n Qo+0.05 

"■■^'^-0.05 


n qo + 0.06 




T2(days) 


537.05 ± 24.00 


f\'?0 1 c;+22. 78 


519 77+29-19 
01»./ '_i3.46 


530.3 ± 27.2 




14.40 ± 1.68 


14 34+1.07 


1 4 21 "'"'^■^■^ 


14.1 ± 1.6 


62 


0.19 ±0.12 


n 1O + 0.09 


o.i4l^;iJ 


0.176 ± 0.174 


■W2 {degrees) 


265.62 ± 45.42 


967 c;q+28.16 


265.28l«;^« 


279 ± 38 


X2 


0.64 ±0.14 


64+°- 1" 
^•"^-0.09 


63+°-^* 




s(ms~^) 


5.47 ±0.69 




^■^°-0.70 




TpiiBJD) 


2453946.95 ± 13.93 


2453946.951^°;°° 


2453950.63l^3;gg 


2453950 ± 10.4 


Ml sin i 


0.87 ±0.06 


r, 07+0. 04 


r, 07 + 0. 06 


0.89 ± 0.12 


ai{AU) 


0.63 ±0.00 


U.Cld_o.oO 


63+° °° 
'^■"'^-O.OO 


0.628 ± 0.02 


Tp2{BJD) 


2454408.22 ± 70.22 


2454408.221^°;°° 


2454420.27l^°|^ 


2454420.3 ± 79.3 


M2 sin i 


0.51 ±0.06 


^•°^-0.04 


"•°^-0.06 


0.504 ± 0.075 


a2iAU) 


1.23 ±0.04 


1 r,q + 0.03 

-^■^•^-0.02 


1 91+0.07 
^■^^-0.02 


1.224 ±0.081 



www.exoplanet.eu, the objective being the reanalysis of the data using a sin- 
gle method and the estimation of corresponding orbital parameters and their 
uncertainties. The results show that while the orbital periods from EXOFIT 
agrees closely with the published ones, the eccentricities show significant vari- 
ation from the published results. This fact indicates that orbital eccentricity 
is not accurately constrained and this degeneracy should be considered more 
carefully. 

We plan to improve to the efficiency of EXOFIT by considering new pa- 
rameterization for the problem as well as faster sampling techniques. Bayesian 
model selection will also be considered for the future versions of EXOFIT. We 
intend to extend our analysis to more planets and planetary data from tran- 
sit photometry to provide a comprehensive Bayesian analysis of the statistical 
properties of orbital parameters of extra-solar planets. 
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Figure 2. Figure on the left plots log(orbital period) of extra-solar planets, 
with results obtained by EXOFIT on the horizontal axis and published ones 
on the vertical axis. Figure on the right shows a comparison of orbital ec- 
centricities with results obtained from EXOFIT on the horizontal axis and 
published ones on the vertical axis. 



References 

Balan, S. T., & Lahav O., 2008, MNRAS, 394, 1936 

Cochran, W. D., Endl, M., Wittenmyer, R. A., Bean, J. L., 2007, ApJ, 665, 1407 
Feroz F., Hobson, M. P., 2008, MNRAS, 384, 449 
Ford E. B, 2005, AJ, 129, 1706 
Gregory, P. C. , 2005, ApJ, 631, 1198 

Jones, H. R. A., Butler, R. P., Tinney, C. G., Marcy, G. W., Carter, B. D, Penny, A. 

J., McCarthy, C, & Bailey J., 2006, MNRAS, 369, 249 
Lewis, A., Bridle, S., 2002, Phys.Rev.D, 66, 103511 
Lomb, N. R. , 1976, Ap&SS, 39, 447 
Scargle, J. D. , 1982, ApJ, 263, 835 



